Motor driven microtubule shape fluctuations - force from within the lattice 
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We develop a general theory of microtubule (MT) deformations by molecular motors generating 
internal force doublets within the MT lattice. We describe two basic internal excitations, the S 
and V shape, and compare them with experimental observations from literature. We explain the 
special role of tubulin vacancies and the dramatic deformation amplifying effect observed for katanin 
acting at positions of defects. Experimentally observed shapes are used to determine the ratio of 
MT shear and stretch moduli (~ 6 x 10 -5 ) and to estimate the forces induced in the MT lattice 
by katanin (10's of pN). For many motors acting on a single MT we derive expressions for the 
end-to-end distance reduction and provide criteria for dominance of this new effect over thermal 
fluctuations. We conclude that molecular motors if acting cooperatively can "animate" MTs from 
within the lattice and induce slack even without cross-bridging to other structures, a scenario very 
much reminiscent of the motor driven axoneme. 

PACS numbers: 87.15.-v 87.16.Ka 87.16.Nn 



Microtubules are the stiffest cytoskelletal component 
and constitute the main routes for motor mediated in- 
tracellular cargo transport in higher organisms [if. Un- 
derstanding their physical properties is at the heart of 
many biological problems from cellular mechanics to in- 
formation and material trafficking in the cell. Since 
the discovery of their high elastic anisotropy [2] it be- 
came increasingly clear that MTs are mechanically more 
complex than other semiflexible biofilamcnts. The high 
anisotropy has been impressively confirmed by thermal 
fluctuation analysis of beads attached to MTs of different 
lengths Q. The emerging picture of the MT is that of an 
anisotropic fiber reinforced material 0, E| with the tubu- 
lin protofilaments (PF) acting as strong fibers weakly 
linked with easily shearable lateral bonds. Remarkably 
this type of design is also found in higher structures like 
axonemes (constituting the backbone of flagella and cilia) 
where relatively inextensible MTs are held together with 
highly stretchable nexin connections [lj]. This remark- 
able structural self-similarity of the two nested structures 
(MT and axoneme) indicates further analogies in the way 
they respond to external and internal forces. We explore 
here important consequences of MT geometry and elas- 
tic properties and show that motors acting on the MT 
surface can generate internal lattice strains sufficient to 
induce observable lateral and longitudinal deformations 
of the MT backbone. 

In the following we describe a twist-free MT of length 
L consisting of N identical PFa with constant dis- 
tance a and a circular cross section, Fig 1. Each PF, 
parametrized by the MT backbone arc length s has a 
position dependent displacement u k (s) from its equilib- 
rium position, k — l,...N. The backbone shape is de- 
scribed by a curvature vector k(s) = 4-t(s) with t (s) 
the bundle center line tangent. The elastic properties 
of the MT are characterized by a PF bending stiffness 
B = ^na 4 Y and compressional modulus K c — ja 2 Y 



withF w 0.1-1.5GPa being the PF Young's mod- 

ulus. Additionally there are shear elastic forces restoring 
the longitudinal displacement between the PFs governed 
by a very soft elastic shear modulus K s w 10~ 3 — IMPa 
[2, [3] ■ The elastic energy is given by 

1 k=N f^/ 2 

£mt-V h / {Bk 2 + K c u' 2 + K S A 2 ) ds (1) 

<5 — /£ - i J-L/2 

with the first term being the bending energy, the second 
the PF compression and the third describing the relative 
shear energy between the neighboring PFs. The shear 
displacement Ak is related to the difference of PF dis- 
placements Uk — Ufe-i and a curvature induced additional 
displacement via 

A fc (s) = u k (s) - u fe _i (s)+ K (s) ■ Ar k ds (2) 

J-L/2 

With Ar'k =f f k- r'k-i and r k = Rmt (cos sin 2**) 
the vector pointing from the MT center to the k-th PF, 
cf. Fig 1. Equations [TJ2] are 3-D analogues of the previ- 
ously proposed stretchable railway-track Q or wormlikc- 
bundlc [5] model for the case of a hollow circular bundle. 
While in general all the N + 3 fields , i.e. the 3 compo- 
nents of k(s) and the PF displacements {u k (s)} fe=1 N 
enter the eqs [T][2] in the limit of small MT deviations 
from a straight line the problem can be drastically sim- 
plified. We first expand the tangent t w (6 x ,0 y ,l) and 

k ks U)' x ,9 y ,0\ in terms of two angular projections 9 X 

and 9 y of t in x and y direction respectively. Exploit- 
ing the circular geometry of the PF arrangement and the 
Fourier representation Uk (s) = u q (s) e « over k we 
quickly realize that only the longest wavelength mode 
iii (s) couples to overall MT backbone shape given by 
the curvature k. This leads to total energy decoupling 
Emt — B\ 1T + Ejtf T + E y MT into a shape- independent 
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FIG. 1: (Color online) The basic geometry of motors induc- 
ing internal force doublets along the MT backbone: between 
two PFs (S-let) and along the same PF (V-let). The red and 
blue ovals represent two coupled motors or two motor sub- 
units (legs) of the same motor. Tubulin lattice vacancies at 
the motor position strongly amplify the MT backbone defor- 
mation (V 3ap -let). 



component E MT and two shape dependent contributions 
(in x and y direction) given by: 
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With A, (s) = a {0, (s) - 9t {-L/2)) - Ui(s) , a = 
\Afk\ the inter-protofilament distance , i = x, y and 
U (s) = {U x ,U y ) = (Rexui,Imxui) with x = 1 - 
e -2-Ki/N an( j renorma ii zec i constants B = NB, K c = 
NK C / (4 -4 cos (2tt/JV)) and K s = NK S . Visually the 
new variable U (s) is a x-y vector at each MT -crosscut 
and can be interpreted as the (vectorial) mean over rel- 
ative PF displacements of neighboring PFs. With this 
enormous simplification at hand we can consider now ba- 
sic motor induced MT excitations (Fig 1). There are two 
elementary configurations in which motors can induce in- 
ternal MT strains: 1) A motor (or a complex of several 
motors) acting between two (not necessarily neighboring) 
PFs and 2) A motor (or a complex of several motors) act- 
ing at two points within the same PF. For reasons that 
will soon become clear we call the excitation 1 S-type or 
simply an "S-let" and excitation 2 we call an V-type ex- 
citation or "V-let". Both excitations are " internal" in the 
sense that there is no net torque or force on the system 
motor+MT similarly to the case of a beating flagellum[lj . 

Elementary internal MT excitations. In the following 
we want to understand the properties of the two basic 



types of excitations from Fig 1 and focus on the S-type 
first. We assume a single motor (or a complex of two 
motors) at position s — Sq bridging between two PF 
with index fej and fc 2 and exerting opposing forces F 
and — F onto them respectively, Fig 1 a) (left) + b). 
The total energy is E tot = E M t + Es-mot with E M t 
given by eqs [TJH and the potential energy of the motor 
Es-mot = -F YlkZkl /-l/2 H s ~ s o) A fe (s) ds . As we 
had for Emt before Es- mo t also decouples into inde- 
pendent modes in the Fourier representation over k and 
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vature independent term and the two shape dependent 



contributions E 



x/y 



Es-mot 



s-mot ( m x an d y direction) given by: 
sq) Ai (s) ds with i = x,y 
and jl being the vector connecting the two attachment 
points of the motor (or motor complex) with components 
Hi = p - ei , Fig 1 b. The equilibrium solution is given 
by the Euler Lagrange equations: 8E\ ot jSUi = and 



with El 



E', 



E l MT and bound- 



J S.mot 

ary conditions 0-(±£/2) = U((±L/2) = (vanishing 
bending and shearing stresses at the ends). A short cal- 
culation leads to A'{ - A~ 2 A l = (1 + a) &£-8 (s - s ) 

with the shear decay length A = (A,T 2 + A^ 2 ) 

Here the two important length scales A c = y K c j K s 



and Xb = yB/a 2 K s have the physical meaning of a 
pure compression- / pure bending- induced shear screen- 
ing length respectively, with their squared ratio a = 
(Ac/As) 2 ~ 35 (N = 13 PF). The remaining equations 
lead to conservation laws (a0i + alii)' = which com- 
bined with the equation for Aj give for the simplest sym- 
metric case so = the tangent angles (up to an arbitrary 
constant) 



Oi (s) = 



cosh(|s|/A-L/2A) 
sinh (£/2A) 



(4) 



With = $f = p ■ ei|4. The resulting MT backbone 
curvature has a jump at s = and attains its maximal 
modulus there. The resulting MT shape is planar (con- 
tained in the plane spanned by jl and t at any position) 
and S-shape-like with initial and final angle coinciding 
Oi (—L/2) = Oi (+L/2) which explains our nomenclature 
"S-type excitation" or "S-let". Interestingly the length 
scale A over which Oi (s) declines allows us to indepen- 
dently estimate the ratio of stretch and shear moduli 
from the observation of S-let deformations coming from 
katanin action^, Y/K s « 647T" 1 (A/a) 2 w 6 x 10 5 for 
A « 1/xm (cf, Fig 2b) and a « 6nm. This value is close to 
the result obtained by Pampaloni et al. Q (Y/ K s « 10 6 ). 
Further the maximal deflection angle of max ~ 38° 
(= 0.66) measured in Fig 2a gives via |4] an estimate 
for the involved motor forces F w 2B# max /(A |/x|) = 
§na 4 Y/{X \n\) = 20-250pA (for \fi\ = a « 6nm and the 
range of values for Y from literature @, Q)- This indi- 
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FIG. 2: a: The shapes of V and S-type excitations as given 
by Eq. [4] and their experimental observations for katanin op- 
erating on MTs (b,c). b: Adapted and edge enhanced image 
from Q. c: Adapted from [(J, black bars highlight the MT de- 
formation. The upper and lower images show the MT before 
and after kink generation respectively. 



cates that many katanin motors might act cooperatively 
to generate the observed shape change in Fig 2b. Another 
interesting possibility is that katanin might be different 
from dynein/kinesin by generating only small contractile 
displacement " powerstrokes" < lnm with a more effi- 
cient chemical-mechanical ATP-energy conversion lead- 
ing to larger contractile forces F > l^ksT/nra = 60 pN. 

The second fundamental internal MT excitation ap- 
pears when a motor (or motor complex) acts along 
a single PF with index k compressing or stretching 
it. The motor energy in this case can be written 
as E V j no t = F (uk (so + H / 2 ) ~ u k (so - \fi\ /2)) « 
F \/j,\ J 6 (s — sq) u' k (s) ds where /i| is the size of the 
motor step. Like in the previous case it is sufficient 
to keep the energy contribution of the mode q = 1 
as the others decouple from each other and from the 
curvature term. Along very similar line of derivation 
as in the S-type excitation case we obtain a solution 
which is planar and contained in the plane spanned by 
the vector r k and t. The resulting tangent angles for 
an excitation in the middle of the MT (sq — 0) are 

given by eqSwith = ^ («) = ft (3 • 
which is now s dependent and changes sign at 
s = 0. Here (j x ,j y ) = sin (n/N) (1 - cos (2tt/N))~ 1 
(sin(27r(A; + l/2)/JV),-cos(27r(A; + l/2)/iV)). The re- 
sulting shape, that we call a "V-let", is a smooth V- 
shaped planar kink in the MT backbone with continuous 
curvature which relaxes on the length-scale A. While su- 



perficially similar the S-let and V-let solutions are physi- 
cally very different for two reasons. First, the s/ \s\ factor 
in the V-let solution changes the symmetry with respect 
to the S-let which leads to dramatic effects on the end- 
end distance as we see below. Second difference lies in 
the different scaling of the numerical prefactors which in 
the case of a V-let do not contain the screening length 
A and involve additionally a very large reduction factor 
1/ (1 + a) . In practice this suppresses significantly the 
involved deformations (9 w 10 -5 — 1CP 6 ), orders of mag- 
nitude below that of a S-let corresponding to the same 
force (6 rj 10 _1 ). However the situation changes dramat- 
ically if the motor is operating at a position of a vacancy 
in the tubulin lattice, cf. Fig 1. In this case the mo- 
tor is not hindered by the large rigidity of the short PF 
portion that the motor acts on. Formally there is no re- 
quirement of continuity for u& (s) of the involved PF k 
at the position of the gap. For such a combined defect 
+ motor excitation which we call V 9ap -let, the motor en- 
ergy is given by E VgaP ^ not = F (u k (s + 0) - u k (s - 0)) 
and its U x / V dependent component becomes E\ 



Vn a P,mot 

= Fji (Ui (s + 0) — Ui (s — 0)) . After a short calcula- 
tion in direct analogy with the previous cases we obtain 
the same form as in eq |4] but with = §>Y ? ( s ) — 
ft (e*j • Arfe) X ~^ F .While its functional form and symme- 
try coincide with that of the V-let, the prefactor of a 
V gap -let solution is more similar in magnitude to that 
of a S-let. It is intuitive to think of a V gap - let as a 
V-let with an effectively renormalized motor force F = 
2a _1 A (1 + a) F. The distribution of angles of observed 
V-shapes measured by Davies et al [(| of max ~ 15° - 
30°(cf.Fig 2c) would suggest for the defect free case (V- 
let) a very large required force of F — 6 x 10 4 — 6 x 10 5 pN . 
However , the same estimate for the V 9ap -\et case gives 
much more moderate 10 — 100piV (few tens of motors), 
showing the prime importance of lattice defects in the 
case of V type excitations. Interestingly Davies et al. Q 
also suggested a crucial role of lattice defects based on the 
pattern and kinetics of MT decomposition by katanin. 

Statistical mechanics of multiple excitations. It is 
particularly interesting to understand the collective 
contribution of a large number of internal MT ex- 
citations acting at random positions Sj and orien- 
tations fj,j in addition to the MT thermal fluctu- 
ations. The motor energy for S- and V-lets be- 
comes now respectively E l Smot — — J2j t 1 )^ ( s j) an(1 

troduce the vector 



= — FJ2j I A*] iWl (sj). It is convenient to in- 



t 1 s/v,j 



such that 



Ms, 



and 



(lyj = \fi\jj (i = x,y). For a fixed (but arbitrary) 
distribution of motors the partition functions Z% = 
J DUiD6i exp (- (E\ [T + E l mot ) /k B T) and the correla- 
tion functions (0i,q6i, P ) (q,p ^ 0) can be obtained. Here 
(...) denotes the average over the random motore distri- 
bution and (...) is the average over the thermal noise. 
\6%,q9i,p} decomposes into the sum of a thermal contri- 
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bution (9i iq 9i tP } T — 2}cbTL 1 G{q)5 p ^ q with the propa- 
gator G(q) = (j3q 2 + ^ and a motor con- 
tribution (Oi, q t ,p) s/V = ^s/v{p)^s/v{q)C S /v{q,P) ■ 
For S- and V-let case we have respectively (p) — 
j0$E]G(p) and * v (p) = j0tk) G ^ 
with the motor position and orientation correla- 
tor C S /v (q, P) = ('Ef,it J -s/v. J t 1 s/v.i C0S ( < l s 3) C0S (P s i))- 
Cg/y (q,p) is easily computed for the simplest choice 
of a uniform motor position and orientation distribu- 



tion: P(si) = l/L , 
with an random angle 



fJ>s/v,j = |Ms/v| (cost 



'3> 



sm< 



;) 



with a probability distribution 



P (4>j) = 1/2tt. In this case we obtain the length reduc- 
tion (Ai) /L « i E 9 ((^) + (^L)) which naturally 
decomposes in a sum of a thermal fluctuation term and 
a motor term (Az) = (Az) T + (Az) In the relevant 
limiting case L/A> 1 we obtain for the thermal part 



(H 



G/; 




fc R T 



a 2 i^° 



(5) 



with the large and small scale persistence lengths 
given by I™ = (j3 + a 2 k)j /k B T and 1° = 

(a/(l + a)) 3 B/ksT. Interestingly the term a 2 K 3 
can be formally understood as an intrinsic self-tension 
straightening the MT at small scales. Similar formulas 
appear in different geometries for the railway-track Q and 
wormlikc-bundlc model[5j. The motor dependent length 
reduction for the S-, V- and V gap -\et excitations with line 
density p is given by: 



cs- 



pF 2 p 



2, ,2 



(Az), 



= c v 



L ^a 4 K 2 X 
pF 2 p 2 L (Az) 
a 2 K 2 X 4 ' I 



(6) 



V ga V pF fl L 
= CV„ap (7) 

a*K 2 X 2 



a 2 (l + a)"^/16 « 6 x 10~ 2 , c v 
1 '* 1 A ~ 4 and cysap = 



with cs = 

= 0.18a 2 (1 + a)~* w 1.3 x 10 
0.73a 2 (1 + a)~ 2 « 0.7. Remarkably the S- and V/V 9ap - 
lets show different scaling. In particular (Az) v , vgap / L 

grows with L (in analogy to the first term in the ther- 
mal contribution [5|) while (Az) s /I/ stays length inde- 
pendent. The physical reason for this difference becomes 
obvious from Fig 2, as the relative slack (Az) /L induced 
by a single S-let scales with X/L, while for an y/V gap - 
let it is essentially length independent. For longer MTs 
this effect leads to a strong dominance of V 9Qp -lets over 
S-lets (Az) Vgap / (Az) s ~ (L/X) > 1. Although hav- 
ing the same L scaling the minute prefactor of defect 
free V-lets renders their contribution relatively insignif- 



long MTs (L — 100/im), underlining the importance 
of lattice vacancies transforming a V-let into a V gap - 
let. Another interesting observation is that in all three 



(Al) 



S/V/Vs^p 



oc pF~. Fixing the number of mo- 
tors N mot but regrouping them into N mot /M clusters of 
size M we have p — » M~ 1 p, F — > MF and therefore 



(A! 



oc M, i.e. the slack grows linearly with 



the cluster size. This indicates that cooperativity (po- 
sitional correlation) of motor action can lead to strong 
enhancement of the slack length. 

From Eqsj5][7]we can derive criteria for the dominance 
of motor slack over the thermal slack. For instance us- 
ing the values estimated from Fig 2 b,c for katanin for 
elastic constants from [3] (F = 20pN, X = 1/xm) and 
p = 8nm, L = 20pm we obtain (Az) jh — p/p c with 
p c S « 0.25nm _1 , p c ,vs^v ~ 1.2 x 10 -3 7im . For large 
enough motor densities the katanin action easily domi- 
nates over the thermal slack (Az) /L « 6 x 10 -4 . 

Being evolutionary specialized for MT deformation and 
degradation katanin is likely to be among the strongest 
slack generating motors. We suspect however that classi- 
cal motors like dynein and kinesin might cause less pro- 
nounced but observable effects as well. While dynein 
is known to walk between several PFs, kinesin is very 
strictly following a single one[l[. Our theory suggest that 
dynein should induce moving S-lets, yet with quickly fluc- 
tuating signs which would diminish the effect consider- 
ably. A battery of many kinesins, however, walking over a 
MT region with many tubulin vacancies, would give rise 
to spatially stationary V 9ap -\ets blinking between "on" 
and "off" states. The theoretical and experimental ex- 
ploration of these issues is an interesting future direction. 

The authors acknowledge fruitful discussions with E. 
Frey, C. Heussinger, M. Bathe, O. Campas, J.F. Joanny 
and P.C. Nelson. I.M.K. acknowledges support by the 
Max-Planck Society. 



[1] 



[2] 
[3] 
[4] 

[5] 
[6] 

[7] 



icant (Az)„/ (Az), 



l X- A L 



10" 



even for very 



J. Howard, Mechanics of Motor Proteins and the Cy- 

toskeleton. Sinauer Press (2001); L.A. Amos and W.G. 

Amos, Molecules of the Cytoskeletion, Guilford (1991). 

A. Kis et al. Phys. Rev. Lett. 89: 248101 (2002) 

F. Pampaloni et al. PNAS 103: 10248 (2006) 

R. Everaers, R. Bundschuh, and K. Kremer, Europhys. 

Lett 29, 263 , (1995) 

C. Heussinger, M. Bathe and E. Frey, |cond-mat/070 2097 
L.J. Davis, D.J. Odde, S.M. Block, and S.P Gross, Bio- 
phys. J. 82, 29162927 (2002) 

J.J. Hartman n et. al, Cell, Vol. 93, 277287; Movie at 
[http://valelab.ucsf.e du /images / mov-rhomtsvkat .mov 
with kind permission by R. Vale 



